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ABSTRACT 

The entropy distribution of the intracluster gas reflects both accretion history of the gas and 
processes of feedback which provide a further non-gravitational energy besides the potential 
one. In this work, we study the profiles and the scaling properties of the gas entropy in 24 hot 
(fcTgas > 6 keV) galaxy clusters observed with Chandra in the redshift range 0.14-0.82 and 
showing different states of relaxation. We recover the gas density, temperature and entropy 
profiles in a non-parametric way, just relying on the assumption of a spherically symmetric 
emission in the deprojection of the best-fit results of the spatially resolved X-ray spectral anal- 
ysis. Adding the hydrostatic equilibrium hypothesis, radial profiles are also obtained from the 
deprojection of the surface brightness, allowing to verify whether the hydrostatic equilibrium 
is a tenable hypothesis by comparison with the spectral measurements. We confirm that this is 
the case on scales larger than 100 kpc and discuss the deviations observed in few non-cooling 
core clusters in the inner regions. We show that the entropy profiles are remarkably similar 
outside the core and can be described by simple power-laws with slope of 1.0 — 1.2. We mea- 
sure an entropy level at 0.1 i?200 of 100 — SOOkeVcm^ and a central plateau which spans 
a wide range of value a few— 200 keVcm^) depending on the state of relaxation of the 
source. The entropy values resolved at given fraction of the virial radius are proportional to 
the gas temperature in these hot systems and appear larger at higher redshift once they are 
compared to the local estimates. To characterize the energetic of the central regions, we com- 
pare the radial behaviour of the temperature of the gas with the temperature of the dark matter 
2dm by estimating the excess of energy /S.E = 3/2 k{Tgas ~ Tum)- We point out that AE 
ranges from « in typical cooling-core clusters to few keV within 100 kpc in non-cooling 
core systems. We also measure a significant correlation between the total iron mass and the 
entropy outside the cooling region, whereas in the inner regions they anti-correlate strongly. 
We find that none of the current models in literature on the extra-gravitational energy is able 
to justify alone the evidences we obtained on the entropy, metallicity and gasH-dark matter 
temperature profiles. 

Key words: galaxies: cluster: general - X-ray: galaxies - intergalactic medium - cosmology: 
observations. 



1 INTRODUCTION 

Clusters of galaxies are the biggest virialized structures in the uni- 
verse, that form at a relatively late time (z < 2). In the hierar- 
chical scenario, the cosmic structures evolve hierarchically from 
the primordial density fluctuations, that are amplified and then col- 
lapse and merge to form larger systems under the action of grav- 
ity. The cosmic baryons fall into the gravitational potential of the 
cluster dark matter halo formed in this way, while the collapse 
and subsequent shocks heat up the intra-cluster medium (ICM) 
to the virial temperature (0.5 ^ T < 10 keV). In the adia- 
batic scenario, the gravity is the only responsible of the physi- 
cal properties of the galaxy clusters, once they are rescaled with 
respect to their masses and epoch of formation jEvrard & Henrvl 
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X-ray properties of galaxy clusters show, however, some 
deviations form this scenario, breaking up the self-s imilarity pre- 
dicted by the adiabatic model (see recent work in Amaud et al.l 
( l2005h : Ponahue et aj] ( l2006h : fonman et al. (20030 and refer- 
ence therein). In parti cular, in the last ye a rs, the studies about X - 
ray scaling relations jEttori et all (|2004): IV ikhlinin et al] l l2005h : 
iKotov & VikhliniiJ l l2005h : iMaughan et al.. ( 2006)) and observa- 
tions of the entropy profiles ^Ponman et al.l (119991) : IPonman et al.l 
(i2003i) ) in groups and clusters of galaxies, and the analysis of 
simulated sourc es with an extra non-gravitational energy injection 
l lBorganil ( l2004l) ) have suggested that we have to account for further 
non-gravitational feedback beyond the gravitational energy. 

The gas entropy records the thermodynamic history of the 
ICM as the product of both gravitational and non-gra vitati onal pro- 
cesses, shaping its observed structure accordingly jVoitl (2005)). 
The measurements of the gas entropy at 0.1 R200 (hereafter 
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5*0.1) showed that the observed value of S is higher than the 
expected one from the adiabatic scenario ( |Ponmanetai]|l999l : 
iLlovd -Davies et al ]|2000*), where S should scale simply with the 
mean temperature of the virialized systems. Instead, an excess in 
the entropy, with respect to the prediction of the adiabatic model, 
is observed in the inner regions of groups and poor clusters at 
some fraction of 7?2oo- This excess sets a minimum value of the 
entropy, labeled as entropy "floor" or "ramp", associated to the am- 
bient gas. The presence of this minimum level of entropy calls for 
some energetic mechanism, not referable to the gravity only, that 
falls into three main classes: preheating, where the gas collapsing 
into the dark matter potential well is preheated by some sources , 
before clusters were assembled at an early epoch ('Kaiser 199l|; 
iBalogh et alJl999l : lTozzi & Normanl200ll;lBorgani et al..200a : lo- 
cal heating by, e.g ., AGN activity, star formation or supernovae 
teialek et al."'200r; 'Srigh enti & MathewdllOO^ : ISabul et al.|[20o3 : 
[Sorgani et al. 2002); cooling, which seems to be able to remove 
low-entropy gas in the centre of th e clusters, producing a similar 
effect to non-gravita tional heating l lBrvanll2000l : iMuanwong et al.l 
l2002l ; lBorganl2004l) . 

In the present work, we aim to confront the model of pre- 
heating, feedback and cooling with the observed properties of the 
gas and of the dark matter in X-ray luminous galaxy clusters, by 
putting constraints on the sources of non-gravitational heating. To 
d o that, we have considered the sample of 24 clusters presented 
in lMorandi et al.l ( l2007h . hereafter Paper I. To quantify the excess 
of energy stored in the ICM with respect to the amount available 
from the gravitational potential, we compare the gas and dark mat- 
ter temperature profile and measure the energy feedback as a func- 
tion of the radial distance. Moreover, the clusters in our sample 
span a wide range of redshift (0.14 ^ 2 ^ 0.82) and have dif- 
ferent state of relaxation. We can thus investigate the dependence 
of the extra-gravitational energy feedback on the cosmic time of 
differently evolved structures. 

This paper is constructed in this way: in Sect.|2] we describe 
the X-ray data reduction and analysis; in Sect. [3] we present our 
results about the entropy distribution in our hot (fcTgas > 6 keV) 
clusters, studying its relation with the gas temperature and metallic- 
ity, its radial profile and how it relates to excess of energy measured 
by comparing gas and dark matter temperatures. We discuss our re- 
sults in Sect.|4]and summarize our findings in Sect. [5] 

Hereafter we have assumed a flat ACDM cosmology, with 
matter density parameter Qom ~ 0.3, cosmological constant 
density parameter SIa ~ 0.7, and Hubble constant Hq = 
70 km/s/Mpc. Unless otherwise stated, we estimated the errors 
at the 68.3 per cent confidence level. 



2 THE DATASET AND THE ANALYSIS 

In Paper I, we describe our dataset and the analysis applied to study 
their X-ray and Sunyaev-Zel'dovich properties. Here, we remind 
the main characteristics of the sample and of the X-ray analysis 
adopted to recover the radial distribution of the ICM entropy inves- 
tigated in the present work. 

We consider 24 galaxy clusters in the redshift range 0.14— 
0.82, emission-weighted temperature between 6 and 12 keV and 
X-ray bolometric luminosity L > 10^^ erg/s, with exposures avail- 
able in Chandra archive. Assuming a spherically symmetric emis- 
sion, the electron density and temperature profiles are obtained by 
deprojecting both the surface brightness profile put in hydrostatic 
equilibrium with a functional form of the dark matter (DM) pro- 



file and the best-fit results obtained in the spatially-resolved X-ray 
spectral analysis by fitting a single thermal component. In partic- 
ular, from the surface brightness profile resolved in a number of 
radial bins between 24 and 239, we obtain directly from the geo- 
metrical deprojection the electron density rij in each j-lh spheri- 
cal shell. The deprojected gas temperature, Tj{q, Pq), is obtained 
by integration of the hydrostatic equilibrium equation once a func- 
tional form of the dark matter density profile, p — p(r, q), is 
assumed, where q =(scale radius, concentration parameter) and 
the gas pressure Po at the X-ray boundary -Rspcc are free param- 
eters. To parameterize the cluster mass distribution, we have con- 
sidered two DM mod els: the universal density profile proposed by 
Navarr o et al. (1997) (hereafter NFW) and the one suggested by 
Rasia et al. C2004) (hereafter RTM). In this study, we adopt the 
RTM model. Our results are not affected if a NFW functional form 
is used. To constrain the 3 free parameters (q, Po), we define a grid 
of values and proceed with a minimization of the merit function 
that compares the observed temperature profile with the projec- 
tion of r, (q,Po) by applying the spectroscopic-like temperature 
definition jMazzotta et al.l2004h . The best-fit values of (q, Po) are 
the ones corresponding to the minimum x^, Xmin- The associated 
errors are estimated at the 68.3 per cent confidence level and are 
computed by looking to the regions in the parameter space where 
Ax^ = ~ Xmin is smaller than a given threshold, fixed accord- 
ing to the number of degrees of freedom (e.g., Ay^ = 1,2.3,3.5 3 
for 1, 2 and 3 d.o.f., respectively; see Paper I and IPress et aljl992h . 
The value of (q, Po) and the related errors are quoted in Paper I. 

Furthermore, we deproject the best-fit results of the X-ray 
spectral analysis, spatially resolved in a lower number of bins (be- 
tween 4 and 10) than the surface brightness profile as requested 
from the higher counts statistic needed to constrain adequately the 
measurements of the temperature. However, in each fc-th shell, the 
electron density Uk and temperature Tk are then recovered without 
any assumption of the hydrostatic equilibrium and provide a direct 
verification of the validity of this assumption once they are com- 
pared to the measures of Uj and Tj described above. The spectral 
deprojection of the observed projected temperature Tproj has been 
performed in a set of n annuli selected to collect at least 2000 net 
counts by inverting the following equation: 

Tproi = {V#{TknlTr)) I {V#{nlTr)), (D 

where the operator # indicates the matrix product (rows by 
columns), V is the effective volume described in Appendix of Paper 
I, and Q = 0.75 using the spectroscopic-like temperature definition 
dMazzotta etal]|2004\ 

2.1 Cooling core and Non-cooling core clusters 

In the following analysis, we divide our sample in three categories, 
depending on the strength of the central cooling-core (see TablelTJ: 

• Strong cooling core (SCC) clusters are the 6 objects in which 
the central cooling time is significantly less than the age of the uni- 
verse at the cluster redshift (tcooi/iagc.z < 0.1). They show very 
low central temperature (~ 2 keV) and strong spike of luminosity 
in the brightness profile, and a very pronounced drop of the tem- 
perature near the boundary of the observation, about a factor two 
compared with the peak of the temperature. The temperature profile 
is very regular, suggesting a relaxed dynamical state. 

• Intermediate cooling core (ICC) clusters have a central cool- 
ing time with values 0.1 < icooi/iagc.z < 1- The six objects in 
our sample show a less prominent spike of brightness than SCC 
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Figure 1. The radial profiles for the projected temperature Tspcc (r). normalized using the cooling-core corrected temperature Tew, for the gas and dark matter 
density are shown for all objects of our sample in the left, central and right panels, respectively. The dashed lines refer to the intermediate cooling core clusters 
(ICC), the solid to the strong cooling clusters (SCC), and the dot-dashed to the non-cooling core clusters. 




Figure 2. Comparison of the entropy and temperature profiles in the internal regions for, from the left to the right, ZW3146 (SCC), AI914 (NCC) and A2218 
(NCC). The points represent each of the measure of Sj in the j-th spherical shell by applying the analysis described in Sect.|2] while the gray region refers to 
the 1 — tr error band. The points with errorbars (triangles) are the measure of S by applying the spectral analysis (see Sect. 12. 2t with (without) applying the 
spectroscopic-like temperature definition of (Mazzotta et al. 2004) . The two triangles on the x-axis refer to the <5 = 0.1 and 5 = 0.3 (see Sect.[3). 



clusters and a mild drop of the temperature in the cooling region 
1/2 Te„). 

• The Non-cooling core (NCC) sources (12 objects in our sam- 
ple) have central cooling time higher than tagc.z and do not present 
any evidence of the central drop in the temperature profile. Both 
the temperature profile and surface brightness map are less regular 
than the ones observed in CC systems, showing hints of substruc- 
tures and merging activity. 



The gas and DM density profiles (right panel of Figure [TJ 
have similar slopes over the entire radial range in the SCC clusters, 
whereas less self-similarity is present in the ICC and especially in 
the NCC clusters: the gas density profile is here flatter than the poM 
one, supporting the scenario in which the ICM has been affected by 
some form of non-gravitational energy. We discuss the physical in- 
terpretation of these observational results in Section [331 



The high level of relaxation of the SCC sources is also con- 
firmed by the study of the polytrophic index 7 0, that has values 
near 1 with a very low scatter for the SCC sources, whereas is more 
scattered in NCC sources at r > 0.57?25oo: 7scc = 1-01 ± 0.09, 
71CC = 1-06 ± 0.12, 7NCC ^ 1.08 ± 0.32. Within 0.3i?2oo, 
we measure 7scc ~ 0.66 ± 0.07, 71CC = 0.97 ± 0.05 and 
7NCC ~ 1-29 ± 0.50, with a clear increase as a function of the 
morphological type and a very high scatter for the NCC sources. 

2.2 On the gas entropy profile 

In the present paper, we have extended the above analysis by esti- 
mating the entropy profile in each cluster by using (i) the gas pres- 
sure Pj and density Uj profile in the equation Sj = PjjrC'J^ and 

^ 7 is calculated as = d log {Tf^)/d log (n^) + 1 by linear fit in the 
log (n/j ) — log (Tfe ) plane by considering the spectral deprojected density 
Wfc and temperature Tj. described in Sect. [2] 
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Table 1. Properties of the sample analyzed. For each object, the name, the redshift z, the emission-weighted temperature Tew, the ratio tcool/tagc, a flag for 
the presence of a strong cooling core, an intermediate one or absence of a cooling core (labeled SCC, ICC and NCC. respectively) are indicated. The columns 
6 and 7 refer to the best fit parameters So.\ and a for the eq.|3]by setting Sq = 0. The last four columns refer to the best fit parameters So, So.i, a and the 
total with the number of degree of freedom (d.o.f.) for the eq.[3] For the redshift and for each of the best fit parameters we report the average value at the 
bottom, by grouping the sources into SCC, ICC and NCC clusters, while the relative errors refer to dispersion of the average values. 
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(ii) the deprojected spectral results Sh ~ Tk/n 
the entropy profiles are obtained by error propagation of the uncer- 
tainties on the single quantity and/or best-fit parameters. We note 
that the dependence of Sj over Po, the gas pressure value at the 
X-ray spectral boundary, can be checked by comparing it with the 
entropy measured once Po is fixed to the value measured in the 
spectral analysis: we find a totally negligible variation at 0.1 i?200 
and a change <. 5 per cent at 0.3 i?200- 

In Figure |2l we present a comparison of the entropy and tem- 
perature profiles recovered with the two methods in the inner re- 
gions of three representative cases, ZW3146 (SCC), A1914 (NCC) 
and A2218 (NCC). We obtain good agreement between the entropy 



measurements in ZW3146 and A1914, the former being an example 
of a typical CC source where the profile decreases moving inward, 
whereas the latter shows the most evident case of flattening, with 
hints of inversion, in the central entropy values. Given the good 
agreement between Sk and Sj even in the internal regions, we be- 
lieve that this inversion is not due to our approach, but it is real in 
A1914 (similar behaviour is found in A773). We note that, if we use 
the entropy recovered by using the proper cooling function in eq. 
[T] instead of the functional T'°', this inversion is less pronounced 
(see Fig.|2ll for Sk- 
in A2218 (NCC), we observe a marginal disagreement be- 
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Sk appears flatter. We draw similar conclusion for A370, A520, 
A2163, and RXJ2228+2037. Nevertheless we observe that for the 
latter sources the low spatial resolution of Sk in the central regions 
(> 150 — 200 kpc) does not allow to sample properly Sj on scales 
of < 50 — 150 kpc, where the inversion occurs. 

We note that the larger deviations between Tj and Tk are ob- 
served in NCC clusters within 100 kpc, where we expect higher 
relative contribution from non-thermal effects due to, e.g., merg- 
ing activity. Therefore, even though the most prominent substruc- 
tures identified in the cluster images were masked, implying that 
we have reduced their effects in the temperature reconstruction un- 
der the hydrostatic equilibrium equation, the sampled gas might 
be still subjected to ongoing merging processes. The higher value 
of Sj compared to Sk in the cluster centre is likely due to a very 
flat density profile that induces a higher temperature value (once 
the hydrostatic equilibrium equation is applied) than the spectral 
deprojected temperature. Indeed, unresolved mergers could lead to 
this very flat density profile (they are clearly visible in A520 and 
A2163), if the gas at i? < 50 — 100 kpc is not wholly relaxed and 
in hydrostatic equilibrium. For the other NCC sources, that do not 
show clearly ongoing merging processes, nevertheless we noted a 
disturbed morphology, as indicated, for example, from the fact that 
the centroid of symmetry does not coincide with the peak of bright- 
ness. 

On larger scales, that involve larger cluster volumes, local de- 
viations from the hydrostatic equilibrium are washed out even in 
the most unrelaxed objects, making tenable the hypothesis upon 
which Sj is obtained. This is also confirmed from (i) the agreement 
between Sj and Sk {Tj and Tk) measured in these sources, (ii) 
the re sults of hydrodynamical numerical simulations jRasia et all 
l200d ). and (iii) the analysis presented in Paper 1 (Sect. 4.1.1), where 
we show how the relation between Aftot and the mass-weighted 
temperature for our sample is in agreement with the results com- 
ing from simulations including feedback and radiative processes, 
supporting our overall mass and temperature Tj reconstruction. 

In the following analysis, we evaluate the entropy at 0.1 and 
0.3 -R200, i-e- at radii well beyond the region where the central in- 
version of Sir) is observed in few NCC objects. Given that, and the 
good agreement on larger scales between the reconstructed profiles, 
we define Sir) = Sj hereafter to fully exploit the spatial resolution 
available in estimating the entropy radial profile. 



3 ENTROPY AND TEMPERATURE DISTRIBUTION 

We examine the S — T relation at fixed overdensities, comparing 
our results with the ones available in literature for nearby systems. 
We investigate, then, the radial entropy profile, studying the be- 
haviour of its gradient and its dependence upon the state of relax- 
ation of the system. Finally, we implement an analysis of the tem- 
perature of the ICM and of the DM to quantify the excess of energy 
associated to the gas and its radial distribution. 

3.1 The entropy-temperature relation 

We have determined the entropy-temperature relation at different 
fraction 5 of the virial radius -R200 {5 = 0.1 and 5 = 0.3). We have 
fitted a power-law model of the form: 

eI'^Ss = aTow,7 , (2) 

where = [J7m(1 + ^)-' + (1 - - nA)(l + z)'' + Q.k] '''^ 
and Tew,? is the total cool-core corrected (by masking the cen- 



tral r — 100 kpc region) emission-weighted temperature in units 
of 7 keV (see Paper 1). The fit has been performed by adopting 
the BCES (Bivariate Correlated Errors and intrinsic Scatter) Y\X 
method ( Akritas & Bershady 1996) (see Paper I for further details 
on this approach). We quote our best-fit results in Table|2]and show 
the distribution of the entropy values at different fractions of -R200 
in Figure [3] We note that SCC clusters show higher normaliza- 
tion (~ 440 and 1400 keV cm^ at 5 = 0.1 and 0.3i?200, respec- 
tively) than ICC and NCC objects, with a larger deviation in the 
inner region {5 = 0.1) which can be explained invoking differ- 
ent relaxation states of the clusters as discussed in Sect. 13.41 The 
best-fit slopes, within the error-bar at la, are in agreement with 
the self-similar prediction (A = 1) and steeper than the slopes 
of A ^ 0.5 — 0.6 observed in l ocal samples of galaxy groups 
and clusters jPiffaretti et al]|2005l : iPratt et al.ll200^ : iPonman et al.l 
I1999II2OO3I) . For comparison, we present in Table |2] also the nor- 
malizations me asured by fixin g A = 1 (self-similar expectation) 
and A ^ 0.65 ( Ponman et al. 2003) and plot in Figur e [3] the best- 
fit results obtained by .Pratt et al., C2006.) and Piffaretti etal.l j2005h 
from their analys es of relaxed groups and clusters at low redshift. 
IPratt e t al.' (20oS) measure A = 0.49 ± 0.15 (a = 271 ± 20) and 
A = 0.64 ± 0.11(a = 990 ± 55) for 5 = 0.1 and 5 = 0.3, re- 
spectivel y, with a clear departu re from the self-similar expectation 
(A = D. jpiffaretti etal] ( l2005h at 5 = 0.1 measure a = 255 ± 71 
by fixing A = 0.65. Once these results are compared with what we 
measure in our sample of very massive systems, we observe that our 
normalizations are on average higher by 20-60 per cent, with slopes 
that are steeper and closer to th e self-similar predicti on than the val- 
ues measured locally (see also lPonman et al. Il l2003h ). This result is 
in agreement with the fact that we are measuring the entropy dis- 
tribution in massive clusters with cool-core corrected temperatures 
in the range 6 — 12 keV. These systems are definitely less affected 
from extra-gravitational, feedback processes that, on the contrary, 
are so relevant in groups and low-mass cl usters representing the 
bul k in the sample of obj ects studied in, e.g.. |Piffaretti et al" I( l2005h 
and lPonman et al.l ( |2003|) . Moreover, by parameterizing the evolu- 
tion in redshift using a (1 -I- z)^ dependence (see Paper I for further 
details on this approach), we did not observe any hints of evolution 
of the entropy-temperature relation within our sample. A very sig- 
nificant positive evolution (B ~ 2 ± 0. 1 but with a reduced of 5 ; 
see Fig.[3]l is instead measured in the relation between the entropy 
estimated at O.I-R200 and Tow when our CC (SCC-l-ICC) objects are 
compared to the best-fit local results in Pratt et al. (2006). Although 
the local best-fits refer to objects distributed over a wider range in 
temperature, the systematic larger values measured at higher red- 
shift is noticeable and definitely more evident at O.I-R200 than at 
0.3i?2oo where we measure _B « 1 ± 0.2 with a reduced of 
about 1. 

3.2 Properties of the entropy profiles 

To characterize the gas entropy profile, we follow 'Donah ue et al.l 
[2006) and fit two different models. The first one reproduces the 
radial entropy profile with a power law plus a constant 5*0: 

Sir)=So + So.i(j--^) (3) 

VO.l 7-200 / 

In the second functional form, we set 5*0 ~ 0, modeling the entropy 
profile with a pure power law. The best fit parameters on the radial 
profile are determined by applying the statistic to the eq.[3]over 
the radial entropy profile between 0.1 7?200 and 0.3 R20Q, whereas 
BCES(Y|X) is used when So is fixed to zero and the region within 
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Figure 3. The S — T relation a different fraction 5 of R200- ^ = 0.1i?200 (left panel) and <5 = O.3R200 (right panel). In each panel the filled circles 
represent the strong cooling core sources (SCC), the triangles the intermediate cooling core clusters (ICC), while the stars the non-cooling core cl usters. The 
sohd l ine refers to the best-fit relation obtained whe n considering all clusters of our sample, while the dashed one represents the best-fit obtained bv lPratt et alj 
i2006l) and the dot-dashed bv lPiffaretti et aljj2005l) . 



Table 2. Best fit parameters of the S — T relation by applying the eq. |2] 
The sources ai'e grouped into SCC, SCC-HlCC and all clusters 





<5 = 


= 0.1 


<5 


= 0.3 


cc 


A 


a 


A 


a 


/NCC 




keV cm^ 




keV cm^ 


SCC 


0.76 ±0.24 


468 ± 70 


0.91 ±0.19 


1380 ± 154 




1 (fixed) 


440 ± 10 


1 (fixed) 


1409 ± 146 




0.65 (fixed) 


494 ± 20 


0.65 (fixed) 


1567 ± 155 


SCC-i-ICC 


1.18 ±0.16 


380 ± 28 


1.17 ±0.22 


1202 ± 94 




1 (fixed) 


399 ± 14 


1 (fixed) 


1186 ±85 




0.65 (fixed) 


432 ± 15 


0.65 (fixed) 


1242 ± 89 


all 


1.33 ±0.20 


354 ± 20 


1.41 ±0.26 


1023 ± 77 




1 (fixed) 


415 ± 11 


1 (fixed) 


1140 ±50 




0.65 (fixed) 


460 ± 14 


0.65 (fixed) 


1216 ±59 




0.1 i?200 is excluded from the fit because it is strongly affected by 
the cooling process. The outermost bins of the fit are excluded by 
the fit, being noisy and likely affected by systematic errors due to 
subtraction of the noise in the data reduction (see Paper 1). Our 
best-fit results are quoted in Table[T] 

The entropy profiles show a regular behaviour (see Fig. 
once the quantities are rescaled to the characteristic value ^^2500 



Figure 4. Profiles of S/ K2500 as a function of R/ iJ2500 ■ The dashed line 
represents the intermediate cooling core clusters (ICC), the solid the strong 
cooling clusters (SCC), and the dot-dashed the non-coolin g core cluster s 
(NCC). The thick solid line represents the profile of Voi t et alj fcOOSi) , 
S/K2rM0 = l-62(r/r250o)"'^'^ (see their Fig. 1, where we have renor- 
malized their entropy profile from A = 200 to A = 2500). 
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at the overdensity o f 250(0 for adiabatic clusters (see, e.g., eq. 2 
in IVoit et alj Jiool) ). Profiles of CC clusters are similar down to 
the inner resolved regions, whereas NCC systems show large de- 
viations in the central parts. These profiles are well reproduced by 
the functional form with a power-law plus a constant for which 
we obtain a Xred always less than or of the order of unity, apart 
from A2390. In particular, SCC sources show a very tight range of 
values of the entropy pedestal So (So < ISkeVcm^) in agree- 
ment with the value found by Donahue et al., and a power law 
behaviour which is roughly preserved on the entire range of the 
radial entropy profile, even in the cooling region (see right panel of 
Figure I?}. The average slope determined from the second method 
(a — 1.18 ± 0.11) is very similar t o the theoretical value of 1.1 
predicted by Tozzi & NormanI ( 1200 ih by using analytic models of 
shock dominated spherical collapse. Concerning S'o.i, it shows val- 
ues in the range 270 — 600 keV cm^: if we adopt the definition of 
Sioo in Donahue et al. as the normalization at 100 kpc, we have 
Sioo ~ 90—150 keV cm^, mildly lower than than the range found 
by them (Sioo ~ 90 - 240 keV cm^). 

The ICC clusters show higher and wider range of So, with a 
typical value of ~ 30 keV cm^ . The power-law behaviour is pre- 
served just on large scale, i.e. outside the cooling region. The av- 
erage slope is still in agreement with above theoretical predictions 
(a = 1.07 ± 0.16), but it is a little lower than the value measured 
in SCC clusters. 

In NCC objects, we observe a more scattered radial pro- 
file, which is likely self-similar beyond the central regions (~ 
0.5 i?2500 ~ 200 — 300 kpc). In the inner regions, we notice 
a very high dispersion on the entropy pedestal value (~ 80 — 
400 keV cm^), larger than the values found in the CC clusters. The 
average slope is mildly lower than the one determined in the CC- 
only subsample (a ~ 0.95 ± 0.21). 

We point out that a rises by considering NCC, ICC and SCC 
sources, respectively (a^^'^ = 1.18 ± 0.11, a^'^'^ = 1.07 ± 0.16 
and Q^*^*^ — 0.95 ± 0.21 for the power-law model). As we will 
see in Sect. 13.31 this trend is probably due to the effect of non- 
gravitational sources on large scale, which justify the flatter radial 
behaviour of the entropy profile in NCC clusters. 

It is worth noticing the behaviour of the entropy pedestal So.i 
in the different subsamples: Sq,?'^ = 300.7 ± 110.3 keVcm^; 
SoT = 57.8 ± 50.7keVcm2; S^'i^ = 6.0 ± 5.9keVcm2. The 
trend of the gas density and temperature profile (see central and left 
panel of Figure[TJ can justify the progressively greater value of the 
entropy in the inner regions by considering SCC, ICC and NCC 
clusters, respectively. We observe higher normalization of the en- 
tropy in SCC sources (see Figure|4]and the value of the parameter 
A in Table[T](. This behaviour is due to the fact that the SCC sources 
show steeper density profiles, i.e. at the same fraction of -R200, as 
long as we consider radii greater than 0.1 -R200, the density of the 
SCC sources is lower. Even though the temperature profiles in the 
SCC sources are a bit steeper than in the ICC and NCC objects, the 
overall effect is that the gas entropy tends to be higher in SCC clus- 
ters. We note that the unrelaxed morphology of the NCC sources 
can not account for systematic changes in, e.g., the determination 
of Ra- 

On the evolution with redshift of the best-fit parameters of 
eq. [3] we note that only for Sj^*^ we obtain a marginal evidence 
of negative evolution (Spearman's rank coefficient rs = —0.60 
for 22 d.o.f with probability of null correlation p — 0.28). On the 



^ -R2500 is ~ O.25R200, i e. ^ 400 - 600 kpc. 
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Figure 5. Normalized break radius )'t,rcak/''200 as a function of the red- 
shift. 

contrary, a shows a positive evolution for the ICC clusters: rl'^'^ — 
0.90, with p = 0.37, while for the SCC and NCC sources there is 
not apparent evolution. 

We have calculated the weighted average value of the slopes 
of the best-fit parameters of the local sample of clusters determined 
by IPonahue et al ] l l2006h . so as to compare it with our estimate at 
higher redshift: they measure a = 1.00 ± 0.01 (when So = 0) 
and OL — 1.23 ± 0.01 (by accounting for 5*0), while we obtain 
a = 1.10 ± 0.01 and a = 1.27 ± 0.01, respectively, by apply- 
ing their procedur^B These results, confirmed also including in the 
sample the ICC sources, suggest that entropy profiles in nearby CC 
systems are slightly flatter than in CC clusters at higher redshift, 
providing some marginal hints on the evolutionary trends present 
in the entropy distribution. 

In Figure [5] we plot the break radius rbrcak present 
in the entropy profile S{r), i.e. the radius where 5*0 ~ 
5*0.1 (rbreak/(0.1 r2oo))" in eq.[3] as a function of the redshift. We 
found the following average values for rbrcak/'"20o: rf£2k/r2oo = 
0.005 ± 0.004, rf,?£k/^200 = 0.050 ± 0.045 and rZ^,^Jr2oo = 
0.172 ± 0.082. The NCC sources show value of rbrcak/'"200 def- 
initely higher than the CC clusters, defining the scale where the 
non-gravitational energy breaks the self-similarity (rbj.^2i ~ 0.1 — 
0.4r2oo ~ 200 — 600 kpc). We do not observe significant evolution 
for r-broak, except for the CC objects (Spearman's rank coefficient 
of rs — 0.89 , probability of null correlation p — 0.019). 

3.3 Gas and dark matter temperature profiles 

In this section, we define a temperature associated to the dark 
matter component following the m ethod presented in llkebe et al.l 
1 I2OO4I) : [Hansen & Piffarettr ( l2007l) . We define the temperature of 
the dark matter halo, Tdm, as: 

/cTdm = i (cr^ + 20-9) finip (4) 

where /i is the mean molecular weight of the ICM, rrip is the proton 
mass, ae and (Jr are the 1-dimensional tangential and radial veloc- 

^ The errors refer to the average value. 
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ity dispersions of the dark matter. The radial velocity dispersions 
has been obtained by solving the Jeans equation: 

GM{(x) _2 /^ rflnpDM(q) , ding,' , 

-ir- = (, d\nR + -dhJi + 2/5(q) ) - (5) 

where a velocity anisotropy parameter is defined, /3(q) = 1 — 
ag/a^. N-body simulations for a variety o f cosmologies show s 
that P has roughly an universal radial profile jCole & Lacevlll99^ . 
which is given by the following relation: 

wherer„ = r/r2oo(q), and/3™ ^ 0.3-0.5 JCarlberg et aljl997h . 
The dark matter profile is estimated as pdm = Ptot — fJ-mpTig^, 
where ptot and Wgas has been determined from the analysis in Pa- 
per I. We will compare the dark matter temperature to the gas tem- 
perature, Tgaa, recovered by applying the hydrostatic equilibrium 
equation Q. 

We solve eq. [5] for /3m ~ {0, 0.4}, corresponding to the 
case of isotropy of the DM and to the central value of the above- 
mentioned range, respectively, to recover a'^ and therefore Tom- 
As boundary condition in eq. [5] we assume Tdm equal to Tgas at 
Rspcc- We have checked that uncertainties on the DM temperature 
assumption at -Rspcc are almost negligible on the DM temperature 
profile in the inner and central regions (R < -R2500 ~ 0.25i?20o), 
being 7?spcc ~ 0.3 — 0.5ii2oo, making our results up to -R2500 
reliable and not affected from the assumed value at the boundary. 
The errors are estimated by looking to the regions of the parameter 
space that satisfy the condition ^ Xmin < 2.3 after the analy- 
sis described in Paper I. Examples of the gas and DM temperature 
profiles for SCC, ICC and NCC objects are shown in Fig.|6] 

Because only the baryonic component is expected to be prone 
to non-gravitational energy effects in galaxy clusters, the difference 
between Tgas and Tdm, AkT, is a powerful tool to trace the ther- 
mal history of the ICM. We show in Fig. [6] how AkT varies as 
a function of the radius. The NCC clusters show a clear trend of 
AkT, with values always greater than zero: AkT « 1 — 2 keV out- 
side the central region 200 — 400 kpc) and it is a few keV in the 
inner region. Near the cluster observed boundary, DM anisotropies 
might make Tdm roughly in agreement with Tgas, even though 
large statistical errors are present and our boundary condition holds. 
A similar trend is observed in ICC clusters, where a less significant 
disagreement between Tgas and Tdm is however observed. In SCC 
clusters, on the contrary, Tgas is well in agreement with Tdm, es- 
pecially in the inner and central regions. 

We notice here that the strong negative evolution measured 
in the scaling relations between yn and the X-ray/SZ quantities 
presented in Paper I, where yn is the integrated Compton parame- 
ter over a fixed angular distance, supports the observed radial be- 
haviour of AkT. Indeed, the measured SZ effect within a fixed an- 
gular size samples larger physical region at higher redshifts. This 
indicates that the effect of non-gravitational processes is relatively 
more pronounced if the SZ flux is measured within smaller physical 
radii, indicating the physical scale over which the non-gravitational 
processes are more relevant. When we perform, instead, the same 
analysis integrating the Compton parameter within a physical ra- 
dius (as done with the quantity i/a), we observe definitely lower 
negative evolution. 

We have also estimated the global excess of energy AC/2500 
defined in this way: 

Following the notation in Sect. |2] 

'^gas — ^7 and Tgas — Tj- 



rR2soo 3 

AC/2500 = / -AfcT(r)ngaB(r)47rr^ dr (7) 
Jo 2 

We find that AC/2500 > 10''^erg in NCC sources (corresponding 
to about 15-20 per cent of the total thermal energy), that is a factor 
between 4 and 10 higher than the measured excess in SCC clusters. 

We refer to Sect. |4] for a discussion of the observational evi- 
dence presented in these two last sections. 

3.4 Relations between gas entropy and metallicity 

The ICM iron mass is a key observable to constrain the cumu- 
lative past star formation history in galaxy clusters. Its relations 
with other observables such as the cluster optical light, total cluster 
mass, stellar mass and gas entropy, together with its redshift evo- 
lution, allow to study the enrichment processes. Moreover, while 
the production of metals is linked to processes of star formation, 
its radial profile is determined by different physical processes, such 
ram-pressure stripping, galactic winds powered by supemovae and 
AGN activity, merger mechanism (Gnedin 1998). 

Following the work of De Gran di et al.l ( 12063) on local clus- 
ters, we present measured iron abundances in the ICM, their evolu- 
tion with the redshift z and their corre lation with the entropy. We 
adopt the solar abundance ratios from I Anders & Grevessd im^ 
with Zq = Fe/H = 4.68 x 10"^ by number. 

We first have determined measures of projected metallicity 
profiles Zpa = Zpeir) = npe/nn, (in units of Zq, that is the 
solar abundance of iron), where npe and nu are the iron and hy- 
drogen densities (by number) respectively. Notice Zp^"-' has been 
integrated up to i?spcc without masking the cooling region, to com- 
pare our results with the literature available. In our sample, we find 
hints of possible negative evolution with redshift, with Spearman's 
rank coefficient of rs — —0.12 for 22 d.o.f ( probability of null cor - 
relation p — 0.59), in rough agreement with lSalestra et al] ( l2007h . 
whose sample covers a wider range of z. After the deprojection of 
the spectral results (see Paper I), we have calculated the iron mass 
enclosed within a sphere of radius R by integrating the iron mass 
density, ppe, over the cluster volume. The total iron mass in solar 
units can be then written as: 

Z 

Mf,{< R) ^ 4nAp,mHYT- / Zp^r) nn{r) r^dr, (8) 
Al© Jo 

where Apc is the atomic weight of iron and mn is the atomic unit 
mass. To integrate the observed profiles at any radius, we have lin- 
early interpolated the metallicity mass profiles within overdensities 
Rs = 0.05 J?2oo, Rs = 0.1i?2oo and A = 2500, which roughly 
correspond to 100 - 150, 200 - 300 and 400 - 600 kpc for our 
sample, respectively. 

No significant evolution with z of Mpc(< i?250o) (probabil- 
ity of null correlation p = 0.94) and Mpe(< J^o.os) (p = 0.52) is 
measured. Instead, we observe a strong segregation between SSC, 
ICC and NCC sources, with SCC clusters that tend to have higher 
metallicity mass by a factor of ~ 2 within 7?2500 and by an order 
of magnitude within R0.05, which roughly corresponds to the cool- 
ing region. The iron mass excess associated with cool core regions 
could be entirely produced by the brightest cluster galaxy (BCG), 
which is always found at the centre of cool core clusters, via SN- or 
AGN- induced winds (De Grandi et al. 2004). Moreover, we con- 
firm the existence of a correlation between Mpc.s and Ss- In Fig.|7] 
we present the correlation between the Mpc,o.o5 as a function of the 
entropy pedestal So (see Sect. |3.2t : we can see an anti-correlation 
between the two quantities, as expected in a picture where the cool- 
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Figure 6. Temperature profiles of the gas (solid line) and of the dark matter (dashed and dot-dashed line for fSm = and /3m = 0.4, respectively). The error 
bands are represented by the gray shaded region for the gas, and hatched region for the DM for the case where /3m = 0. The clusters are A1835 (SCC), A2261 
(ICC), and A2218 (NCC), from the left to the right. 
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Figure 7. Congelation between Mpo,o.05 and the entropy pedestal So (Isft), Mfo,o.i and 5o.i (right). 



ing is the likely predominant physical process in the cooling re- 
gion. Enrichment from recent Supernovae type la in the cD galaxies 
can explain the central metal abundance e xcess observ ed in cooling 
core clusters (De Grandi et al. 2004; Bohringer et al.i2004) . On the 
contrary, outside the cooling region, we observe again nearly self- 
similar relation between Mpc,o.i and So.i (see Fig.|7) as we have 
seen in Sect. l3.1l for the 5*^ — Tew relation. 

Assuming a synthesi zed iron mass per SNIa event mia of 0.74 
Mq jNomotoetaUliggTh and an energy output of 10 erg, we es- 
timate that 1 — 3 X 10^° SNIa events in the region inside i?2500 
are required to produce MFe,2500 ^ 0.8 x lO^'^M© observed in 
NCC clusters. This number of SNe corresponds to a global energy 
output of 1 — 3 X 10^^ erg over the entire lifetime of the cluster, 
that is lower by a factor 2-4 than the excess of energy Ai725oo esti- 



mated in Sect.[33](AC/'^*^'^ ~ 10''^ erg), suggesting the action of 
other sources of non-gravitational energy to fully account for this 
observed excess. 



4 DISCUSSION 

The main results emerging from our study of the entropy profiles 
in hot (fcTgas > 6 keV) galaxy clusters at 2 > 0.1 are that these 
profiles, although similar in the outskirts where they behave as a 
power-law with slope 1.0 — 1.2, are remarkably discrepant in the 
central regions, with SCC objects that show a power-law behaviour 
down to the innermost spatially resolved regions and NCC clusters 
having profiles that flatten to a constant value at r < 0.3i?25oo 
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(Fig. |4ll. Accordingly, thie comparison between gas and dark mat- 
ter temperature profiles (Fig. [T} reveals that SCC clusters do not 
present any significant energy excess at any radius, whereas ICC 
and, more dramatically, NCC objects show Ai5 = 3/2AkT larger 
than few keV in the cooling region and above. Note that the situ- 
ation near the boundary of the sources is unclear, because the sta- 
tistical errors are very large, the effect of possible anisotropies in 
the DM are there more prominent and we are assuming some con- 
straints on the dark matter temperature at iispec- This excess of 
energy with respect to the 'gravitational energy floor' associated 
to the DM temperature profile can be interpreted as an indication 
of the presence of some form of non-gravitational energy that can 
constrain the mechanisms affecting the ICM thermal history and 
the observed entropy profiles. Indeed, in agreement with AE ~ 
over the entire radial range, we observe that SCC clusters have very 
low entropy pedestal values So of few keV cm^, while the higher 
and more scattered values of So measured in ICC and, particularly, 
in NCC systems can be justified by an injection of energy AE 
of 1-10 keV, that, distributed over scales < 100 — 300 kpc, ex- 
plains also their flatter entropy profiles. The regular behaviour of 
the entropy profiles outside 0.1 f?200 is also in agreement with the 
fact that AE is low at these radii, where we have to consider the 
limitations of our analysis near iispec as mentioned above. This 
scenario is also supported from our results on the S — T rela- 
tion, where we observe an higher normalization, more significant 
in the inner regions (5 — 0.1; see Table[3} of the SCC subsam- 
ple with respect to ICC and NCC sources. We note hints of larger 
entropy values at higher redshift when our measurements in CC 
clusters are compared to the best-fit results obtained in nearby sam- 
ples, with a more significant deviations observed at 0.1 R200 than 
at 0.3-R200, suggesting that cores in our CC objects are not yet 
well defined from the radiative processes. Moreover, the observed 
mild differences in the slopes of the entropy profile, with a that be- 
comes slightly higher by considering NCC, ICC, and SCC sources, 
respectively (a^^^ = 1.18 ± 0.11, q'°*^ = 1.07 ± 0.16 and 
^NCC _ Q gg _^ Q 21 for the powerlaw model) can be explained 
by looking at the temperature and density profiles (Figure[Tll, which 
are a bit flatter for non-cooling core systems: this trend can be justi- 
fied with small energy excess (AE ~ 1-2 keV) at large scale in the 
NCC objects compared to the NCC and SCC ones. The radial be- 
haviour of AE(r) is also confirmed by the analysis made in Paper 
I, where we noticed a strong negative evolution in the j/n— X-ray 
and ynSZ scaling relations (see Sect. l3.3t . 

All our systems are the products of the hierarchical scenario, 
how is suggested from the similar behaviour of the gas temperature, 
density, entropy and dark matter profiles in the regions above the 
cores. On the contrary, the cooling region characterizes SCC, ICC, 
and NCC systems. In particular, continuous interplay between cool- 
ing and some form of (pre-)heating can explain the variety of the 
properties observed, with SCC dominated from the cooling phase 
and, on the other end, NCC still subjected to some effects of heat- 
ing. 

Theoretical models must predict the magnitude of the ob- 
served AE{r), and the impact of the non-gravitational processes 
associated to this excess in the central regions. These models 
fall into three main classes: preheating, where the gas collapsing 
into the dark matter potential well is preheated by some mech- 
anism , before clusters we r e assembled at an early epoch |Kaisei 
199l |: iBalogh et all 1 19991 : iTozzi & Normal I2001I : iBorgani et al" 
20051 ): local heating by AGN activity, star formation or supemovae 



low-entropy gas in the centre of th e clusters, producing a similar 
effect to non-gravita tional heating ( lBrvanll2000l : iMuanwong et al.l 
I2OO2I : lBorgiuill2004h . Hereafter, we review the main characteris- 
tics of these models and discuss how they are consistent with our 
observational constraints. 



(iBialek et alJl200ll:lBrighenti & Mathewi 2006: Ba bul et al.|[20o3 : 
iBorgani et al 2002h : cooling, which seems to be able to remove 



4.1 Preheating models 

Models of pre-heating, where a constant energy input is injected 
either prior of the cluster collapse (0.1 — 0.3 k eV per particle, e.g. 
iNavarro et al] ( ll995h : lTbzzi & Normaiil 1 200 ih) or after the cluste r 
formation (1 — 3 k e V per particle (e.g. Metzler & Evrard ('1994'); 
ILoewensteinl ( l2000h : IWue t al. (2000); Bower etal. (2001), could 
justify, only partially in NCC objects, the observed magnitude of 
Ai?(r), but not its ra dial behaviour. Nevertheless, as pointed out by 
IBorgani et al.l ( l2005h by studying hydrodynamical simulated clus- 
ters, there is no possibility to inject a large quantity of energy per 
particle (< 1 keV), unless a large isentropic core is produced in the 
entropy profile, core t hat is not observed in our profiles in agree- 
ment with other work dPonman et id]|2003l : [Pratt & Arnaudl2003l 

2004 

Ponman et ajj j2003h suggest that any raise of the tempera- 
ture and/or decrease of the density in the gas inside the primor- 
dial structures due to preheating can get largely raised by the ac - 
cretion shock. Following the model of lDos Santos & Dorel ( l2002h . 
IPonman et al.l ( |2003|) estimated that a mild raise of the entropy of 
the gas confined to filaments (~ 10 — 100 keV cm^, coiTesponding 
to a temperature of ^ 10~^ keV) can be boosted by the accretion 
shock to the observed value of S (~ 100 — 1000 keV cm^). They 
point out that an interplay between shock and smoothing of the pri- 
mordial gas due to a preheating can justify the observed properties 
of the gas entropy, given the above upper limits on the energy bud- 
get of the preheating and being the slope of the entropy profile close 
to the value predicted from shock heating. 

Preheating prior of the cluster collapse should be a energeti- 
cally favorable mechanism compared to in situ heating to cast fur- 
ther energy into the gas before it is concentrated in the gravitational 
potential well of the DM halo, since less energy is required to in- 
crease the entropy of the gas by a given amount when its density 
is lower as in the filaments. In fact we observe that in the shock 
dominated collapse scenario, a mild injection of energy through 
preheating can greatly amplify the final energy Ea^ of the post- 
shoc ked particles, b eing Enn oc Ein, with _Ein the initial energy. 
Borg ani et al. 

liQ03) 

show that smoothing the accretion pattern by 
preheating in the case of simulations without radiative physics am- 
plifies the entropy generation out to the radius where the accretion 
shock acts. Nevertheless, the effect seems to be substantially re- 
duced when cooling is also taken into account. 

However, our estimates of AE{r) show not a constant pro- 
file but instead a declining one outwards. This behaviour can not 
be explained by any preheating mechanism, either prior or after 
cluster collapse, even though entropy is amplified through subse- 
quent shock heating. Results for X-ray bright nearby objects by 
[Pratt et al. ( 2006) support this conclusions, because their scaled en- 
tropy profiles show increasing scatter in the inner regions, with a 
dispersion (~ 60 per cent) definitely higher than the value fou nd in 
simulations including filamentary preheating (~ 30 per cent, IVoitI 
120051) ). 
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4.2 Heating models 

In principle, heating can amplify the boosting of the entropy out 
to the radius where accretion shocks are taking place, especially in 
low mass system, since they are accreted by smaller subhalos where 
the gas is more smoothed by the extra heating. As pointed out by 
iBorgani et alj dlOOS}) . local heating due to star formation activity 
is not able alone to prevent overcooling a nd to reproduce the pre- 
dicted star formation as low as measured jMuanwong et alj|2002l) 
and the observed entropy profile: maybe we have to appeal to fur- 
ther sources of non-gravitational energy, like AGN, not taken into 
account in such hydrodynamical simulations, or different physical 
mechanisms to distribute the energy inside the ICM. 

The need of this further source of non-gravitational energy is 
also confirmed by the analysis made in Sect. 13.41 where we ob- 
serve that the number of supemovae we require to reproduce the 
observed metallicity is not able to account for all the excess of en- 
ergy AC/2500- 

A gentle, transonic heatin g process, such as t he weak shocks 
detected in the Perseus cluster jpabian et al1l2003l) . can provide a 
framework by which one can explain all the observed properties, 
like the flattening of the entropy profile in the innermost regions (<: 
a few tens of kpc) even of SCC clusters, and the spikes of metallic- 
ity measured in the centre of SCC sources (see Fig.[7}. Weak shocks 
are indeed likely not able to prev e nt me t als' ac cumulation in the in- 
nermost regions. [Ponahue et alj j2005l . [200^ pointed out that the 
central cooling time of the SCC galaxies (~ 10* yr) is consistent 
with the time scale of the activity of radio sources (~ a few 10^ 
yr) at the centre of clusters. Energy casted by the radio jet 10*^ 
erg/s) can then produce the observed flattening of the entropy pro- 
files on scale of a few tens of kpc. Gasdynamical models of jets 
flows proceeding from a central supermassive black hole and en- 
tering surrounding gas may heat the ICM by casting mass and en- 
ergy outwar ds till scale ~ hundreds of kp c, possibly lowering the 
cooling rate jSrighenti & Mathew j|2006l) . Nevertheless the above 
picture does not explain the excess of energy Ai5 in the regions 
outside the core, where we have probably to require some other 
form of heating, like e.g. shocks induced from merging activity. 

Whatever sources of non-gravitational energy we have to ap- 
peal, they must fuel energetically the ICM in such a way to re- 
produce the magnitude and the radial behaviour of A£'(r), casting 
energy (in the order of few keV) and metals preferentially in the 
innermost regions. 



4.3 Cooling models 

Cooling plays a key role to explain the observed excess of energy 
A_E(r). In fact, Ai5(r) is higher in the central regions moving 
from SCC to ICC and to NCC systems. Cooling can easily account 
for this trend through radiative losses of the accumulated thermal 
energy. In the SCC clusters, Tgas is roughly equal to Tdm sug- 
gesting either that a perfect balance between cooling and heating 
is established, permitting the radiative losses of the only amount 
of energy in excess with respect to the one associated to the DM, 
or that heating is episodic and we are observing structures in their 
undisturbed phase. 

In a similar manner, IVoit etaljj2002l) have argued that the en- 
tropy responsible for the break of the self similarity is not a global 
property of the ICM, but rather a property set by radiative cooling: 
they point out that the observed entropy value at the core radii of 
groups and clusters is near to the entropy at which tcooi ~ tnubbic. 

A model to explain the observed features in the entropy profile 



is that proposed bv lVoit&BrvanI toll) . They show how cooling 
and supernovae heating act to eliminate high compressible gas with 
S < Sc, being Sc the cooling threshold, from the X-ray phase. 
Much of the condensation and the feedback is prior of the epoch of 
clusters' formation, balancing these processes reciprocally: more 
the cooling is effective, more the star formation is active with re- 
lease of energy to the ICM and consequent reduction of the cooling 
itself. This picture is likely not wholly adequate, leading to a very 
large isentropic core in the entropy profiles, which are not observed. 



5 SUMMARY AND CONCLUSIONS 

We have presented Chandra observations of the entropy profiles 
and scaling properties of a sample of 24 galaxy clusters spanning 
the redshift range 0.14-0.82 and classified accordingly to their cen- 
tral cooling time in strong (SCC), intermediate (ICC) and non- 
cooling core (NCC) systems. We have performed a spatially re- 
solved spectral analysis and recovered the gas density, temperature 
T and entropy S profiles at high spatial resolution and in a non- 
parametric way. We have shown that those entropy profiles are re- 
markably similar outside the central regions with a typical entropy 
level at . 1 7?200 of 1 00 — 500 ke V cm^ , and have a central entropy 
plateau covering a wide range of values (~ a few— 200 keV cm'^), 
with the highest values associated to NCC objects. The CC clusters 
show larger values of the entropy than the one measured in nearby 
luminous systems, with a more significant deviations observed at 
0.1 -R200 than at 0.3 -R200 with respect to the best-fit results in Pratt 
et al. (2006), suggesting that the core in our CC objects are not yet 
well defined from the cooling processes. 

We have studied the radial behaviour of the temperature of the 
gas (Tgas) and of the dark matter (Tdm). We have found that Tgas 
is always higher than Tdm: for the SCC clusters, the difference of 
temperature AfcT — fcTgas — fcToM is negligible, while it is large 
for the non-cooling core clusters (up to -^10 keV per particle), with 
AkT{r) that declines outwards. 

We conclude that none of the models of (pre-)heating and 
cooling discussed in literature is able to explain alone the observed 
trends of the entropy profiles and of AE{r). Likely, we have to re- 
quire an interplay of these processes. A scenario with shock dom- 
inated collapse and preheating in the primordial filaments might 
account for most of the extra-gravitational energy, as confirmed by 
the slopes of the entropy profiles near th e theoretical value of 1. 1 
expected in the accretion shock picture (Tozzi & Norman* l200ll) . 
As described in Sect. 14.11 this should be an energetically favor- 
able mechanism compared to the in situ heating, amplifying sig- 
nificantly the final energy Efi^ of the post-shocked particles. On 
the other side, gentle, sub-sonic heating processes, e.g. supported 
from AGN's jets, can account for many of the observed properties, 
but not for excess of energy still measured outside the core (see 
Sect.[331). 

The relaxed environment of SCC clusters is required to en- 
hance the central metal abundance and total iron mass in correspon- 
dence of low entropy regions (see Sect. |3.4| l. There, the cooling is 
so efficient to remove on short time scale the excess of energy per 
particle of the ICM, permitting (i) Tgas to approach the dark matter 
value Tdm and (ii) the other physical parameters used in describing 
the entropy profile, like, e.g.. So, rbroak, to vary. 

Overall, the different observed behaviour of the entropy pro- 
files of SCC, ICC and NCC massive clusters suggest that we are 
observing the end products of the hierarchical model for structure 
formation. They represent different stages of the relative relevance 
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of heating and cooling in regulating the feedback that shapes the 
ICM distribution: galaxy clusters are identified either as NCC ob- 
jects when heating, probably due to a residual merging activity and 
feedback from AGNs triggered from the merger itself, is predom- 
inant, or as sec systems when the radiative losses are energeti- 
cally prominent, being ICC objects an intermediate case between 
the two. 
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